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The precision of the CP phases 2/3 and 2/3 s determined from the mixing-induced CP asymmetries 
in Bd —>• J/ipKs and B s — » J/ip<fi, respectively, is limited by the unknown long-distance contribution 
of a penguin diagram involving up quarks. The penguin contribution is expected to be comparable 
in size to the precision of the LHCb and Belle II experiments and therefore limits the sensitivity of 
the measured quantities to new physics. We analyze the infrared QCD structure of this contribution 
and find that all soft and collinear divergences either cancel between different diagrams or factorize 
into matrix elements of local four-quark operators up to terms suppressed by Aqcd /mg,, where my, 
denotes the J/i/> mass. Our results, which are based on an operator product expansion, allow us to 
calculate the penguin-to-tree ratio P/T in terms of the matrix elements of these operators and to 
constrain the penguin contribution to the phase 2/3 as |A<£d| < 0.68°. The penguin contribution to 
2/3„ is bounded as |A<$}| < 0.97°, |A0s| < 1.22°, and |A<^"| < 0.99° for the case of longitudinal, 
parallel, and perpendicular <j> and J/iji polarizations, respectively. We further place bounds on 
|A0d| for Bd —> ip(2S)Ks and the polarization amplitudes in Bd —> J/t/K*. In our approach it is 
further possible to constrain P/T for decays in which P/T is Cabibbo-unsuppressed and we derive 
upper limits on the penguin contribution to the mixing-induced CP asymmetries in Bd —> J/ipn°, 

Bd —» J/ipp°, B a — » J/t/Ks, and B s —> J/ipK*. For all studied decay modes we also constrain the 
sizes of the direct CP asymmetries. 

PACS numbers: 13.25.Hw 

INTRODUCTION 


The mixing-induced CP asymmetry in the decay Bd —> 
J/ipKs is the key quantity to measure the CP phase of 
the Bd — Bd mixing amplitude. Within the Standard 
Model (SM) this CP asymmetry Aq$C j ^ Ks (t) deter¬ 
mines the angle /3 = arg[— V t bV t * d /(V c bV* d )\ of the unitar- 
ity triangle. The B factories BaBar and Belle had been 
designed to measure A^C J ^ Ks (t) to a high precision to 
probe the Kobayashi-Maskawa (KM) mechanism of CP 
violation. Within the Standard Model, the KM phase 
is the only source of CP violation in weak transitions 
and therefore must correctly describe all CP asymme¬ 
tries measured in weak hadron decays. The measurement 
of (3 at the B factories gave us sufficient confidence that 
the KM mechanism correctly describes CP violation in 
both K and Bd decays and led to the dedication of the 
2008 Nobel Prize in Physics to Makoto Kobayashi and 
Toshihide Maskawa. Today’s focus of flavor physics is 
the search for physics beyond the Standard Model which 
reveals itself in small deviations from the KM picture. In 
generic models of new physics B—B mixing probes new 
physics associated with scales beyond 100 TeV; reducing 
the uncertainties of Standard-Model predictions is there¬ 
fore of utmost importance. B s —> J/ip4> is the analogous 
key mode in the B s —B s system. Since the unitarity of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix essentially 
fixes /3 S = arg \-V// b V ts / {y c bV* s )\ = 1.0° to a very small 
value, CP studies of B s —► J/ifcj) directly probe physics 


beyond the SM. The decay amplitude Af for an b —> ccs 
decay B q —► /, where / is a CP eigenstate consisting of 
a charmonium state and a light meson, can be written as 

Af = X s c Tf + X s u P f (1) 

with A® = V* b V ps , p = u,c , and 

Tf = %U\ c 'Q° + c *Q*+X' CiQi \ Bq) ' (2) 

V 3 

Pf = ^(f\C 1 Q u 1 +C 2 Q^+Y / C :Qj\ B q )- ( 3 ) 

Here, Q\ = s“t m (1 - 75)q ,/3 9 /3 7 M (l - 7 b)b a and Q 2 = 

s a 7 M (l — 75)<7 a <f 3 7 /i (l — 75)6 /3 are the current-current op¬ 
erators. The index j labels the penguin operators Qj 
which involve the CKM elements Af = — A* — A®. While 
the QCD penguin operators Q 3-6 and Q$g are impor¬ 
tant for this paper (see Ref. [l] for their definition), elec- 
troweak penguin operators have negligible effects. The 
time-dependent CP asymmetry A C p J (t) = [T(B q (t) —> 
f)-T(B q (t) -7 f)]/[T(B q {t) -7 f)+T(B q (t) -7 /)] reads 

A B ^f (f ) = S f sin(A M q t) - C f cos(A M q t) 

CP cosh(Ar 9 i/2) + Hai", sinh(Ar g t/2) 

Here A M q and Ar g are the mass and width differ¬ 
ence, respectively, between the mass eigenstates of the 
B q —B q system. We write — 77 / sm(<f> q + A^> q ), where 

CP I/) = r/f\f) and <f> q is the CP phase in the limit 
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Pf = 0. The SM predictions are 4>d = 2/3 and </> s = —2 (3 S . 
To first order in e = \V us V u b/{y cs V c b)\ « 0.02 one has 

Pf 

tan(A</>) ~ 2esin 7 Re—. (5) 

T f 

Comparing Eq. © (with 7 = (69.7 ± 2.8)°) with the 
present experimental world average sin cpd = 0.679±0.020 
| 2 | (meaning an error of 1 . 6 ° for cpd) shows that the pen¬ 
guin contribution already matters now and will certainly 
do so for future measurements at LHCb and Belle II. Tf 
and Pf are non-perturbative multi-scale matrix elements, 
which defy calculations from first principles of QCD. 

For the prediction of the branching ratio B(B d —> 
J/ijfKs) one only needs Tf, which was addressed with 
the method of QCD factorization Q in Ref. Q: in the 
limit of infinite charm and bottom masses Tf can be 
expressed in terms of the J/ip decay constant and the 
Bd —> Ks form factor. The result of Ref. Q underesti¬ 
mates B(Bd J/ipKs ) by a factor of 8 . This failure, 
however, is not surprising, because the corrections to the 
infinite-mass limit are of order Aqcd /(m c a s ) and there¬ 
fore numerically unsuppressed for the actual value of the 
charm mass ; 3| [5|. The standard approach to quantify 
Pf/Tf in Bd —> J/4’Ks uses the approximate SU(3)f 
symmetry of QCD (or its U-spin subgroup) which relates 
the decay of interest to b —> ccd modes like B s —► J/t/jKs 
and Bd — > J/i/m 0 0 , 0 - A drawback of this method is 
our poor knowledge of the quality of the SU(3)f sym¬ 
metry in Bd tS —> J/ipX (with X = Ks, 7 r 0 ,...) decays. 
(Comparisons of branching ratios essentially test SU(3)f 
in Tf only, with little sensitivity to Pf.) Furthermore, the 
b —> ccd control channels have 20 times smaller statistics 
than their b —> ccs counterparts. SU(3)f seemingly fails 
in B s —> J/ip4 > : because the </> meson cannot be closely 
approximated by an SU(3)f eigenstate, but is an equal 
mixture of octet and singlet. 

In this paper, we present a dynamical calculation of 
Pf/Tf which does not assume an approximate SU(3)f 
symmetry. Our results permit, for the first time, the 
prediction of Sf and C/ also for b —> ccd decays. 


OPERATOR PRODUCT EXPANSION 

For definiteness we first specify the discussion to Bd —>■ 
J/i/Ks and return to B s —> J/ip(p and other modes in 
the phenomenology section. For B{B d J/ipKs) we 
only need Tf and can neglect the penguin coefficients. It 
is useful to express T f in terms of the matrix elements of 

Qov = s7 M (i - 75)607^0, 

Qoa = S7 M (1 - 75)^07^750, 

Qsv = S 7 M (1 ^ 7 5 )T a 6c7 M T a c, 

Qsa = S 7 M (1 - 7 5 )T a 6c7 /i 7 5 T a c. (6) 


Then Tf in Eq. © becomes Tf = ^(J/tpKs\Co(Q 0 v - 

Qoa ) + Cs(Qsv ~ QsA)\B d ) with C 0 = C 2 /N c + Ci and 
Cg = 2 C 2 , where N c = 3 is the number of colors. Using 
next-to-leading order (NLO) Wilson coefficients in the 
naive dimensional regularization (NDR) scheme [l|, [sj at 
the scale /r = m^ one finds Co = 0.13 and Cg = 2.2. 
The smallness of Cq is a well-known numerical accident 
entailing that the weak decay produces the (c,c) pair 
almost in a color octet state. We normalize the matrix 
elements (for j = 0, 8) as 

( QjV) = VoVj, ( QjA) = Vo aj (7) 

to the factorized matrix element Vq = (Qov)fact = 
2fj/^m Bd PcmF^ K {ml) = (4.26 ± 0.16) GeV 3 . The 
uncertainty stems from the form factor F/ S ^ K /rn/,) = 
0.586 ± 0.021 0 and the J/\b decay constant fj = 
(0.405 ± 0.005) GeV. m Bd = 5.28 GeV and p cm = 
1.68 GeV are the Bd mass and the three-momentum of 
the J/ip or Ks in the Bd rest frame. Uo,8,ao,8 depend 
on p in such a way that the /z-dependence of Co, Cg can¬ 
cels from physical quantities. When we quote numerical 
values we refer to the choice p = m^. The large- N c 
counting of our (complex) hadronic parameters is vq = 
1+0(1/Nl), vg, ag = 0(1/N C ), and a 0 = 0(l/N/). Nor¬ 
malizing the branching ratio to the experimental value we 
find 

B(B d ->■ J/^Ks) 

B{B d -A J/lpKs)ex p 

[1 ± 0.08] |0.47u 0 + 7.8(vg - a 8 )| 2 . (8) 

Varying the phase of vg — ag between —n and 7 r one finds 
the correct branching ratio for 0.07 < |us—as] < 0.19 if uo 
is set to 1. Thus, there is no mystery with the branching 
ratio and the hadronic parameters obey the hierarchy 
expected from 1/N C counting. The terms involving ao are 
negligible in view of other uncertainties and are omitted 
throughout this paper. 

Pf in Eq. © receives contributions from Q “ 2 and the 
penguin operators Qj, j > 3. The matrix elements of the 
latter can be trivially expressed in terms of the operators 
in Eq. ©• Therefore, this contribution to Pf/Tf only 
depends on vg/vo and ag/vo■ Below we will see that the 
magnitudes of these ratios are under control thanks to 
the 1 /N c hierarchy of vo,Vg, ag and the information from 
B(B d -A J/t/Ks)ex p- By varying the parameters in the 
allowed ranges we can then find the maximal contribution 
of the penguin operators to |A</>|. 

In order to apply the same strategy to 2 we must 
first express the up-quark penguin depicted in Fig. [la] 
in terms of matrix elements of the local operators in 
Eq. ©. In Ref. 0 it is argued that a penguin loop 
flown through by a hard momentum q (in our case 
q 2 ~ = (3.1 GeV) 2 ) can be calculated in perturbation 

theory (“Bander-Soni-Silverman (BSS) mechanism”). In 
Ref. [llj this idea is used to find an estimate of (Q%) 
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FIG. 1: The LO diagram is shown in (a). The soft IR divergence of the diagram (b) factorizes with the 
corresponding diagram of the effective-theory side shown in (c). The diagram (d) is an example of a diagram with a 

collinear IR divergence. In (e) a spectator diagram is given. 


which leads to an upper bound on |A0| which is smaller 
than the values found by SU(3 )f arguments [3|- In 
this paper, we turn the BSS idea into a rigorous field- 
theoretic method by proving an operator product expan¬ 
sion (OPE) 

(J/4’K s \Q]\B d ) = J2 C jtk (J/^Ks\Qk\B d ) + • ■ - (9) 

k 


with k running over k = OF, OA, 8V, 8A. The dots repre¬ 
senting terms suppressed by higher powers of Aqcd / yfq 2 - 
The Wilson coefficients Cj t k = C j°J + (a s (/j,)/(47r+ 
... are calculated in perturbation theory to the de¬ 
sired order in a s (p), with the renormalization scale 
p = 0(m^,mb). A similar OPE has been derived to 
calculate charm-loop effects in the rare semileptonic de¬ 
cays B K^£ + £~ 12]. Since leptons carry no color 
charges, this application involves no four-quark opera¬ 
tors like those in Eqs. © and (0. From Fig. [la] one finds 


c: 


2 = 0 except for C 8G,sv = - 
with the penguin function 

2 a 


( 0 ) 


and <%% = P(q 2 ) 


P(q 2 ) = 


3 47t 


In (4) 

2 

— in — — 

L \d 2 ) 

3 J 


( 10 ) 


Inherent to applications of the OPE as in Ref. [12[ or 


in this paper is the assumption that rescattering effects 
for values of q 2 far above the partonic pair-production 
threshold are correctly described in perturbation theory. 
Eq. 0 captures all hadronic effects in the (u,u) —► (c, c) 
transition only if there is no intrinsic ( u,u ) component 
in the J/ij) wave function (e.g. no J/ip-p° mixing). A 
powerful check of our framework will be the confrontation 
of our predictions for b —> ccd transitions with data. 


PROOF OF FACTORIZATION 

In order to establish Eq. 0 we must prove that the 
coefficients Cy*, are infrared (IR) safe. To this end we an¬ 
alyze i) the soft IR divergences of the two-loop diagrams 


contributing to (Q“), ii) the collinear IR divergences of 
these diagrams, iii) spectator scattering diagrams, and 
iv) higher-order diagrams in which the large momentum 
bypasses the penguin loop (“long distance penguins”). 


i) An example of a diagram with a soft divergence is 
shown in Fig. EB This soft divergence is reproduced 
by the corresponding diagram of the effective-theory side 
(i.e. RHS) of Eq. 0, depicted in Fig. EB so that this 
divergence factorizes with Cj°2 and does not affect . 
All soft divergences are from diagrams in which the addi¬ 
tional gluon connects two external lines and cancel from 
Cj 1 ^ in the same way. 


ii) Collinear divergences occur in diagrams in which 
a gluon is attached to the line with the strange quark, 
which we treat as massless. An example is shown in 
Fig. Ed] If l denotes the loop momentum flowing through 
the gluon propagator and p s is the momentum of the 
external strange quark, the collinear divergence corre¬ 
sponds to the region with l 2 = 0 and l oc p s . We can then 
reduce the problem to the study of one-loop diagrams 
with an external on-shell gluon: If we sum over all possi¬ 
bilities to attach this gluon to one of the lines of the LO 
diagram in Fig. E3 the collinear Ward identity of QCD 
ensures that this sum vanishes when the open Lorentz in¬ 
dex of the gluon line is contracted with P. This feature 
ensures that the collinear divergences of the sum of the 
two-loop diagrams vanish. (For a discussion in the con¬ 
text of QCD factorization see Refs. 001 0].) It equally 
holds for the effective-theory side of the OPE. The can¬ 
cellation of collinear divergences is conceptually identical 
to the situation in typical processes in collider physics; 
it is further known to be much simpler (with fewer dia¬ 
grams to be discussed) if a physical gauge (with only two 
propagating gluon degrees of freedom) is adopted. 


iii) Next we discuss the spectator scattering contribu¬ 
tions: diagrams in which the gluon connects the b or s line 
with the spectator quark line trivially factorize with the 
corresponding diagrams on the effective side. If the gluon 
connects the spectator with the gluon line or a charm or 
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up line, we have to take into account that the squared mo¬ 
mentum in the penguin loop is ( q+l ) 2 instead of q 2 . If the 
gluon is soft, P ~ Aqcd, the expansion of the loop func¬ 
tion P around q 2 reproduces a term which correctly fac¬ 
torizes with up to term suppressed by Aqcd / \fq 2 ■ If 
the gluon is hard-collinear, with virtuality l 2 ~ PchiAqcd, 
where p cm ~ 1.5 GeV is the three-momentum of the Kg 
or J/if) in the Bd rest frame, the situation is more sub¬ 
tle: the LO diagram is suppressed by Agco/Pcm, because 
the momentum of the spectator quark changes from zero 
to 0(p cm ) in the decay, which is penalized by the light- 
cone distribution amplitude (LCDA) of the kaon Q. The 
asymptotic form of the kaon LCDA, <I>(:r) = 6x(l — x), 
where x and 1 — x are the fractions of the kaon momen¬ 
tum carried by the s and d quarks, favors momentum 
configurations in which the kaon momentum is roughly 
equally shared between the two valence quarks. While 
the propagator of the scattered hard-collinear gluon is 
suppressed as ~ l/(AQCDPcm), the suppression of the LO 
diagram is lifted, because the spectator momentum is in 
the region x ~ 1/2 favored by the kaon LCDA. To iden¬ 
tify further suppression factors we first discuss the case 
that the gluon connect a charm line with the spectator: 
counting q 2 ~ m ^ and the energies of s and spectator-d 
quarks as p cm / 2, the penguin loop gives P((q + l) 2 ) ~ 
P{q 2 ) + The non-factorizing piece involving 

the derivative P'{q 2 ) comes with a factor of Pcm/vn^. The 
virtuality of the (anti-)charm propagator is around p cm 
entailing a suppression factor of Aqcd/Pciii- Thus, alto¬ 
gether spectator scattering from the charm lines obeys 
Eq. ([9]) up to terms of order AqcdA*V- Next we dis¬ 
cuss the spectator scattering from the up line, with a 
sample diagram depicted in Fig. [lc] We find that these 
diagrams are power-suppressed by Aqcd /m^. In this 
respect these spectator diagrams differ from the similar 
photon penguins calculated in Ref. fl5| . which involve 
P{q 2 ) for q 2 ~ 0 rather than q 2 ~ rri 2 ,. 

vi) So far we have assumed that the underlying hard 
process is the penguin loop with the hard scale s/q 2 . But 
it may also be possible that the hard momentum trans¬ 
fer to the J/ip occurs through a hard gluon radiated from 
the b or s line, while the penguin loop is a “long-distance 
penguin” governed by soft QCD. Such a situation is ex¬ 
emplified by the diagram in Fig. [lb] with the left gluon 
having virtuality ~ m^,. These diagrams, in which the 
whole weak decay process occurs with small momentum 
transfers, have a suppression factor (Aqcd /x/q 2 ) 3 stem¬ 
ming from the hard gluon propagator and an off-shell b 
quark propagator (or s quark propagator). 

In our power counting in i)-iv) we have treated p cm 
as an intermediate scale between Aqcd and and 
have found no non-factorizable non-perturbative effects 
of order p cm /m While p cm enters two-loop diagrams 
through pb ■ p s ~ rribPcm ., such terms do not come with 
IR divergences and end up in the NLO corrections to the 


coefficients Cy*,. We find that the counting rule for p cm 
is irrelevant, one can reproduce our results above as well 
for the limiting cases p cm ~ Aqcd and p cm ~ yfq 2 ■ In 
particular, higher orders of the OPE do not involve op¬ 
erators with derivatives acting on the s field. The same 
feature was found for B —>• K in the last paper of 

Ref. 0. 

The choice q 2 = m ^ for P{q 2 ) may be altered by 


adding a contribution of order Aqcd to \[q 2 . This 
shuffles a piece proportional to (Aqcd /rn^)P'(rri^) into 
the coefficient of the sub-leading operator S 7 M (1 — 


ls)T a b 


O-ml 


C 7 a ‘T“c, which removes the ambiguity 


associated with the choice of q 2 . At NLO in a s one 
generates non-zero coefficients cj 1 ^ also for j = 1 or 
k = OA, 8A. 

In conclusion the OPE with the minimal set of opera¬ 
tors in Eq. works, the coefficients in Cy*, are IR-safe. 


PHENOMENOLOGY 


The penguin amplitude depends on the Wilson coeffi¬ 
cients as 

P f = Ro (2C 4 + 2C 6 + 2 C 2 C^ 8V + C 8G di% 8V ) v 8 + ... 

( 11 ) 

where the dots represent the terms with vq and a 8 
which have much smaller coefficients. The dependence 
of C^° 8V (calculated from Fig. fTal) on the renormalization 
scheme cancels with the scheme dependence of C 4 + Cq 
in Eq. CD. In the NDR scheme adopted by us these 
penguin coefficients give a larger contribution to Pf than 
the u-penguin loop contained in C, 2 ° 8V . This is not sur¬ 
prising, because the u-penguin loop enters at NLO, while 
C 4 + Cq already contributes at LO. The omission of this 
dominant LO piece explains the smallness of the result 
in Ref. 111. 

For the prediction of Pf/Tf we implement the con¬ 
straint from B{B —> f ) exemplified for / = J/ipKg in 
Eq. @ in the following way: adapting a phase conven¬ 
tion in which Af in Eq. m is real and positive, we can 
determine a 8 in terms of Vo, vq, vg, and the measured 
B(B —> f) |l 6 j, [17]. Then we use this to eliminate a 8 
from Pf/Tf. For example, we find 


Pj/i!>Ks 

Tj/^Ks 


= 0.01 -0.02u o - (0.71 + 0.33i)v 8 (12) 


for the central value of Vo (quoted after Eq. 0) and 
B(Bd —► J/tpKs). We vary v 8 and vq in their allowed 
ranges |t>s| < 1/3 and |vo| = 1 ± 0.15 with the constraint 
that |cX 8 1 < 1/3 must be obeyed. The allowed ranges for 
A(f>, Cf , and AS/ = Sf +rjf sin <p q are almost symmetric 
around zero. We list the upper bounds on their mag¬ 
nitudes for several decay modes in Tabs. |T] and [0 The 
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results include the uncertainties from Vo, the branching 
ratios, CKM parameters [l8j], and higher-order terms in 
our OPE. For the b — > ccd decay modes with Cabibbo- 
unsuppressed Pf/Tf the expansion in Eq. ® has been 
replaced by the exact formula (see e.g. Ref. 1 , 0 )- Our 
bounds are conservative, as the considered ranges for 
vs and as are wide (permitting even sizable cancella¬ 
tions in Eq. ([8])). From Eqs. (0 and (fl2l) one veri¬ 
fies that any additional information on magnitude or 
phase of one of these parameters will substantially re¬ 
duce the ranges quoted in Tabs. Q] and |TTJ Our results 
for Bd —y J/'ipn 0 favor the Belle measurement Cj/^o = 
—0.08 ±0.17, Sj/^no = —0.65 ±0.22 [if} over the BaBar 
result Cj /,^0 = —0.20±0.19, Sj/^o = -1.23±0.21 0. 
(In the absence of penguin pollution Cj/^ n o = 0 and 
Sj /^ro = — sin(2/l) = —0.69 ± 0.02.) In the the case of a 
more precise and non-vanishing measurement of Cj/^o, 
for example, Cj/^o = — 0.10± 0.01, which corresponds 
to the current world average with a ten times smaller 
error, we can also put stronger restrictions on the shift 
of the mixing-induced CP violation \ASj/^°\ < 0.13. A 
measurement of Cj/^q that is consistent with zero, how¬ 
ever, does not improve the bound. This feature occurs 
in all decay modes with Cabibbo-unsuppressed Pf/Tf. 
The measurements of Sf and Cf for the Bd —>■ J/il>p° 
polarization amplitudes [21] comply with the ranges in 
Tab.[U 


CONCLUSIONS 


We have established a factorization formula (to lead¬ 
ing power in Aqcd /m^,) for the penguin contribution to 
the CP- violating coefficients S / and C / in A c p 1 (t) for 
final states / containing charmonium and the related 
shift Atj> q of the corresponding CP phase. As a cru¬ 
cial result the penguin contributions involve the same 
hadronic matrix elements as the tree amplitude. This 
allows us to constrain Pf/Tf, which determines Sf and 
Cf, and to find e.g. |A</> d | < 0.68° for B d —*■ J/'ipKs 
and |A</i^| < 0.99° for B d —> representing bounds 

that were thought to be uncalculable from first prin¬ 
ciples. Novel territory are our predictions for Sf and 
Cf in b —> ccd decays, in which Pf/Tf is Cabibbo- 
unsuppressed. Future experimental probes of these pre¬ 
dictions will constitute a powerful test of our theoretical 
framework, whose key ingredient is an operator product 
expansion for the up-quark penguin loop. There are no 
similar consistency checks for the standard predictions 
of Pf/Tf based on SU(3 )f symmetry, which, moreover, 
cannot be used for B s —»■ J/ipcj). We further remark that 
our results do not depend on any properties of the char¬ 
monium LCDA. 
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